Triple approach to determination of the c-axis penetration depth in BSCCO crystals 
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The c-axis penetration depth \ c in Bi2Sr2CaCu20s+a (BSCCO) single crystals as a function of 
temperature has been determined using three high-frequency techniques, namely: (i) measurements 
of the ac-susceptibility at a frequency of 100 kHz for different sample alignments with respect to the 
ac magnetic field; (ii) measurements of the surface impedance in both superconducting and normal 
states of BSCCO crystals at 9.4 GHz; (iii) measurements of the surface barrier field H.j{T) oc 1/ X C (T) 
at which Josephson vortices penetrate into the sample. Careful analysis of these measurements, 
including both numerical solution of the electrodynamic problem of the magnetic field distribution 
in an anisotropic plate at an arbitrary temperature and influence of defects in the sample, has 
allowed us to estimate A c (0) ~ 50 /jm in BSCCO crystals overdoped with oxygen (T c « 84 K) 
and A c (0) « 150 /xm at the optimal doping level (T c w 90 K). The results obtained by different 
techniques are in reasonable agreement. 
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I. INTRODUCTION 



Study of anisotropy of the magnetic field penetration 
depth as a function of temperature in high-T c supercon- 
ductors (HTS) advance considerably our understanding 
of pairing-mechanism in these materials. It is known (see, 
e.g., Ref.EJ and references therein) that in-plane penetra- 
tion depth A\ ab (T) oc T in the range T < T c /3 in high- 
quality HTS samples at the optimal level of doping, and 
this observation can be interpreted the most simply in 
the ci-wave model of the high-frequency response of HTS. 
Measurements of out-of-plane or c-axis penetration depth 
A C (T) are quoted less frequently than those of X a b(T). 
Most of such data published so far were derived from mi- 
crowave-measurements of the surface impedance of HTS 
crystalstrH There is no consensus about AA C (T) at 
low temperatures. Even in high-quality YBCO crystals, 
which are the most|-^>tudied objects, one can find rboth 
linear, AA C (T) oc TKra, and quadratic dependences!!] in 
the range T < T c /3. In BSCCO materials, the shape 
of AA C (T) depends on the level of oxygen doping: in 
samples with-Kaaximal T c ~ 90 K AA C (T) oc T at low 
temperaturesa'Q; at higher oxygen contents (overdoped 
samples) T c is lower and the linear function AA C (T) 
transforms to a quadratic oneQ. The common feature of 
all microwave experiments is that the change in the ratio 



AA c (T)/A c (0) is smaller than in AA a b(T)/A a b(0) because 
in all HTS A c (0) > A afc (0). Another possibility to deter- 
mine c-axis penetration depth is the measurements of the 
first penetration field of Josephson vortices. In quasi-2D 
systems, their penetration may be impeded by a surface 
barriet- the value of which is inversely proportional to 
A c (T)tLil. The quantitative estimates of A C (T) deduced 
from the surface barrier data were however disputedta. 
Furthermore, A C (T) is._aJso inversely proportional to 
the plasma frequencytJ'til which is usually assigned to 
Josephson plasma resonance frequency modified by the 
field-dependent interlayer phase coherence^ IJ. How- 
ever, this interpretation is still controversialtj. There- 
fore, independent measurements of A C (T) are of inter- 
est. To date, all the above mentioned properties have 
been studied separately. It is the aim of this paper to 
apply together three different techniques to the deter- 
mination of the absolute value of A C (T) in order to ob- 
tain unambiguous results: (i) ac-susceptibility measure- 
ments of BSCCO crystals have allowed us to determine 
the temperature variation AA C (T); ii) cavity perturba- 
tion technique and electrodynamic analysis of the surface 
impedance anisotropy is used to determine both the vari- 
ations and absolute values of A C (T) and A a f,(T); (iii) the 
first penetration field of Josephson vortices is measured 
and shown to be related to A C (T). 
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II. ELECTRODYNAMIC BASIS OF THE 
MEASUREMENTS 



The electrodynamics of layered anisotropic HTS is 
characterized by the components a a b and a c of the con- 
ductivity tensor. In the normal state, ac field pene- 
trates in the direction of the c-axis through the skin 
depth S a b — and in the CuC>2 plane through 

5 C = ^J2/oj[1q(j c . In the superconducting state all pa- 
rameters S ab , 8 C , cr ab = a' ab - icr" b , and a c = a' c - ia" are 
complex. In the temperature range T < T c , if a' <C cr" , 
the field penetration depths are given by the formulas 

y/ '1 / 'oj \lqg" ■ In the close neigh- 
the decay of magnetic field 



\ib = ^/l/w/xocr" b , A c 
borhood of T c , if a' > a" 
in the superconductor is characterized by the functions 
Re ((5 a b) and Re (S c ), which turn to 5 a b and 5 C at T > T c , 
respectively. 

In this paper we present the results of high-frequency 
measurements of anisotropy of two BSCCO single crys- 
tals with various levels of oxygen doping. The first 
sample (#1), characterized by a lower critical temper- 
ature, T c « 84 K (slightly overdoped), has dimensions 
a X b X c ~ 0.8 x 1.8 x 0.03 mm 3 . The second (#2, 
a x b x c ~ 1.5 x 1.5 x 0.1 mm 3 ) is almost optimally 
doped (T c re 90 K). 
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FIG. 1. Curves of the ac-susceptibility of sample #2 versus 
temperature in different orientation with respect to ac mag- 
netic field: H„ || c (full circles); _L c, H^, is parallel to 
the b-edge of the crystal (up triangles); H u _L c, is par- 
allel to the a-edge of the crystal (down triangles) . Left-hand 
inset: transverse (T) orientation, H„ || c, the arrows on the 
surfaces show directions of the screening current. Right-hand 
inset: Longitudinal (L) orientation, H u _L c. 



For determination of both X a b(T) and A C (T) compo- 
nents of the penetration depth we measured the tem- 
perature dependences of the Q-factor and the frequency 
shift Sf of the resonant circuit for different sample align- 
ments with respect to the ac magnetic field H w : in the 
transverse (T) orientation, H w || c, when the screening 
current flows in the a6-plane of the crystal and in the 
longitudinal (L) orientation, H w _L c, with currents run- 
ning in the directions of both CuC>2 planes and the c-axis 
(insets to Fig. 1). In the first case the values of Q(T) 
and Sf(T) are directly connected with in-plane penetra- 
tion death X a b{T) at T < T c and skin depth S ab (T) at 
T > T C EJ. Both lengths are smaller than the thickness of 
the crystal. In the second L-orientation at T < 0.9 T c 
the penetration depth is still smaller than characteristic 
sample dimensions. But at T > 0.9 T c the lengths A c and 
S c are comparable to the width of the crystal. In order to 
analyze our measurements in both superconducting and 
normal states, we used formulae for field distributions in 
an anisotropic long strip (b 3> a, c) in the L-orientation. 
These formulae neglect the effect of ac-faces of the crys- 
tal, if H w is parallel to the 6-edge of the crystal (inset on 
the right of Fig. 1), but take into account the size effect. 
At an arbitrary temperature, the measured quantities are 
expressed in terms of the complex parameter [i intro- 
duced in RefO, which is controlled by the components 
Vab(T) and c c (T) of the conductivity tensor through the 
penetration depths 8 a b and <5 c t3: 
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where the sum is performed over odd integers n > 0. In 
the superconducting state, if a' <C cr", the parameter \x 
is real: 
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III. AC-SUSCEPTIBILITY MEASUREMENTS 

The first approach to determination of the c-axis pen- 
etration depth A C (T) in BSCCO single crystals is based 
on measuring the ac-susceptibility x — x' — *x" at a fre- 
quency of 100 kHz. The imaginary part x" °f X is propor- 
tional to the energy dissipation in the sample and the real 
part x' is proportional to the shielding of magnetic field. 
Ac-susceptibility is characterized by the components of 
the conductivity tensor and, hence, allows to determine 
the penetration depth A. Experimentally it is possible 
if the penetration depth is comparable to the sample di- 
mension. So the ac-susceptibility measurements are usu- 
ally adapted to the superconducting powders whose grain 
size is comparable to A. However, since A c (0) in BSCCO 
single crystals is relatively large, we managed to deter- 
mine A C (T) from the temperature dependences of x'(T). 
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As an example, figure 1 shows the temperature de- 
pendences x'(r)/lx'(0)l hi sample #2 for three differ- 
ent sample alignments with respect to the ac magnetic 
field: in the T-orientation, H w || c, (full circles); in the 
L-orientation, H w _L c, (H w is parallel to the &-edge of 
the crystal, up triangles); in the L-orientation, whose dif- 
ference from the previous configuration is that the sam- 
ple is turned around the c-axis through 90° (down trian- 
gles). Fig. 1 clearly shows that at T < T c x! ab (T) is no- 
tably smaller in the T-orientation than Xab+J?) m the 
L-orientation (the subscripts of x' denote the direction 
of the screening current). The coincidence of x' a b+c(T) 
curves at H u _L c and the small width of the supercon- 
ducting transition at H w || c (AT C < 1 K) indicate that 
the quality of the tested sample #2 is fairly high. 

In the superconducting state at T < 0.9 T c we find that 
A a & "C c and A c <C a. In this case, we derive from Eq. (2) 
a simple relation between the real parts of ji and x : 



2A a f, 
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FIG. 2. Femperature dependences AA C in sample #1 
(squares) and #2 (circles) at T < 0.9 T c . Open symbols plot 
low-frequency measurements, full symbols show microwave 
data. The inset shows low temperature sections of the AA C 
curves in sample #2. 



Figure 2 shows measurements of AA C (T) in sample #1 
(squares) and sample #2 (circles) at T < 0.9 T c . The 
open symbols plot low-frequency measurements obtained 



from x'(T) in accordance with Eq. (3), the full sym- 
bols plot microwave measurements (see the next section) . 
Agreement between these measurements is fairly good, 
but in fitting together experimental data from sample #2 
(upper curve) we had to divide by a factor of 1.8 all 
AA C (T) derived from measurements of ac-susceptibility 
using Eq. (3). The cause of the difference between 
AA C (T) measured-in sample #2 at different frequencies 
is not quite clearta. 

The curves of AA C (T) at T < 0.5 T c plotted in Fig. 2 
are almost linear: AA C (T) cx T. The inset to Fig. 2 shows 
the low-temperature section of the curve of AA C (T) in 
samale #2. Its slope is 0.3 /im/K and equals that from 
Rcf.El. Note that changes in AA C (T) are smaller in the 
oxygen-overdoped sample #1 than in sample #2. 

We also estimated A c (0) on the base of absolute mea- 
surements of the susceptibility Xc(0) an d we obtained 
A c (0) » 70 fim for sample #1 and A c (0) « 210 /im for 
sample #2. 



IV. CAVITY PERTURBATION TECHNIQUE 

The second technique of determination A C (T) is mea- 
suring the difference A(l/Q) (A/) between reciprocal 
Q's (resonant frequency shifts) of a cavity with a sample 
inside and the empty cavity as functions of temperature 
at a frequency / = 9.4 GHz. These parameters are re- 
lated to the surface impedance Z s = R s + iX s compo- 
nents through the geometrical factoc, r s of the sample: 
R s = r s A(l/Q), AX S = -2r s A///li The penetration 
depth is A(T) = X s {T)/lj^ at T < T c . In order to 
determine the absolute value of X s — —2T s Sf/f, where 
Sf is the frequency shift relative to that which would be 
measured for a sample with perfect screening and no pen- 
etration of the microwave fields, one needs the constant 
parameter f = A/(T) - Sf(T). In HTS this constant 
can be derived from microwave measurements in the nor- 
mal state. 

In the T-orientation /o can be derived from the condi- 
tion that the real and imaginary parts of the impedance 
should be equal above T c (normal skin-effect )E3. Given 
fo, we can calculate the conductivity a a b{T) — u' ab — 
i<j" b = iu> fio / (T) at all T. The temperature depen- 
dences of cr' ab (T) (open squares) and cr" h (T) (open cir- 
cles) are shown in Fig. 3a and Fig. 3b for crystals #1 
and #2 respectively. Note that cr' ab (T) in Fig. 3b docs 
not have a broad peak at low temperatures, typical 
for <j' ab {T) in high-quality HTS. The reason for that is 
the rather large value of the residual surface resistance: 
R ICS = | _2 S (T — > 0) w 0.5 mfi in the a6-plane of the sam- 
ple #2c!J, i? ros was four times less in the sample #1. 
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T(K) T(K) 
FIG. 3. Conductivities a' ab (T) (open squares) of two BSCCO single crystals #1 (fig. a) and #2 (fig. b) at 9.4 GHz, extracted 
from the surface impedance measurements. The insets show a'^ b (T) = A^ b (0)/A^ b (T) data (open circles). 

In the L-orientation of a crystal shaped as a long strip 
the quantities A(l/Q) and 5f are expressed in terms 
of the complex function fj, = fx' — from Eq. (1) or 

A(l/Q)-2iSf/f = i^v/V, (4) 
where v and V = 58 cm 3 are the volumes of the sam- 



ple and cavity respectively, 7 = 10.6 is a constant of 
our cavitytl. In the superconducting state at T < 0.9 T c 
the expression for fj, is given by Eq. (3). The curves of 
AX C (T) in sample #1 (full squares) and sample #2 (full 
circles), measured by cavity perturbation technique and 
obtained using Eqs. (3), (4), are shown in Fig. 2. 
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FIG. 4. Temperature dependences of A(l/Q) (open squares) and — 2A/// (open circles) at H u || b of sample #1 (fig. a) 
and sample #2 (fig. b) at T > T c . Solid lines show the functions —25f(T)/f deriving from Eqs. (4) and (1). Left-hand insets: 
A(l/Q) and — 2A/// as functions of temperature in the neighborhood of T c . Right-hand insets: p c (T) (triangles). 



We can estimate A c (0) by comparing of A(l/Q) and 
A/ = 5f — /o measurements taken in the T- and L- 
orientations to numerical calculations by Eqs. (1) and 
(4), which take account of the size effect in the high- 
frequency response of an anisotropic crystal. The pro- 
cedure of comparison for sample #1 and #2 is illus- 



trated by Fig. 4a and Fig. 4b respectively. Unlike the 
case of the T-orientation, the measured temperature de- 
pendence of A(l/Q) in the L-orientation deviates from 
(— 2A///) owing to the size effect. Using the measure- 
ments of cr a i,(T) at T > T c in the T-orientation (Fig. 3), 
alongside the data on A(l/Q) in the L-orientation (open 
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squares in Fig. 4), from Eqs. (1) and (4) we obtain the 
curves of p c {T) = l/a c (T) shown in the right-hand in- 
sets to Fig. 4a and Fig. 4b for crystals #1 and #2. Fur- 
ther, using the functions <J C (T) and a a j,(T), we calculate 
(—2Sf/f) versus temperature for T > T c , which are plot- 
ted by the solid lines in Fig. 4. These lines are approxi- 
mately parallel to the experimental curves of —2 A/// in 
the L-orientation (open circles in Fig. 4). The difference 
— 2(6 f — A/)// yields the additive constant /o- Given f 
and A/(T) measured in the range T < T c , we also obtain 
Sf(T) in the superconducting state in the L-orientation. 
As a result, with due account of \ a b(T) = \J '1 / 'uj HQ(r" b (T) 
(insets to Fig. 3), we derive from Eqs. (4) and (3) A c (0), 
which equals approximately A c (0) 50 /im in sample #1 
and 150 /<m in sample #2. These results are in reason- 
able agreement with our ac-susceptibility measurements 
if we take into consideration the fact that the accuracy 
of A c measurements is rather poor and the error can be 
up to 30%. The temperature dependences of (t"(T) for 
both BSCCO single crystals are shown in Fig. 5. 
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FIG. 5. Conductivities a"(T) of BSCCO single crystals #1 
(left scale) and #2 (right scale) at 9.4 GHz, obtained by com- 
paring the measurements of A(l/Q) and A/ = 5f — fo to 
numerical calculations by Equation (1). 



V. MICROWAVE ABSORPTION 
MEASUREMENTS IN A STATIC MAGNETIC 
FIELD: SURFACE BARRIER AND INFLUENCE 
OF DEFECTS 

An alternative technique for the determination of A c (0) 
is based on the measurements of the magnetic field Hj (T) 
at which Josephson vortices penetrate into the sample. 

In the experiment the microwave absorption measured 
in the L-orientation as a function of the static magnetic 
field (0-30 Oe) parallel to the CuC>2 layers exhibits a no- 
table increase at the field Hj(TjE!l. Figure 6 displays the 
change of microwave dissipation, starting from zero field 




H(Oe) 

FIG. 6. Microwave absorption as a function of the applied 
field H at 3 temperatures, for 2 orientations (0° and 2°) of 
the H with respect to the afe-plane. The onset of dissipation, 
indicated by the arrow, occurs at the penetration field Hj(T). 



(within ±0.1 Oe) measured in sample #1 for various ori- 
entations of the applied field close to the afr-plane (0° and 
2°) and in a low-field range: < H < 25 Oe, at three 
typical temperatures (T = 78 K, 65 K, 50 K). After each 
field sweep, the sample was warmed through T c and then 
cooled again in zero field, in order to avoid any possible 
vortex pinning when studying the penetration starting 
from zero field. The dissipation of Josephson vortices is 
characterized by the fact that it does not depend on the 
angle (Fig. 6), as long as these vortices remain locked. As 
the field increases, an onset in the dissipation occurs at a 
temperature-dependent field Hj(T), which we associate 
to Josephson vortices entering the sample. As in Ref.c3, 
we choose to define Hj(T) as the field value where the 
microwave absorption exceeds the experimental accuracy. 
The field thus determined is plotted in Fig. 7. The error 
bars take into account both the noise and the estimated 
drift of the signal with time. 
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FIG. 
tal #1 
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Plot of Hj(T) (full circles) for BSCCO crys- 
Triangles display the upper bond of H c i(T) using 
A a f)(0) = 2100 A and A c (0) = 10 pm. The temperature varia- 
tions AX a f,(T) and AA C (T) are taken from the present work. 



The magnitude of Hj(T) is too large to be associated 
with the first thepaodynamic critical field H c i(T) for 
Josephson vorticesca: 



H cl (T) 



47rA afc (T)A c (T) 



[lnA ab (T)/d+1.12], (5) 



where <po is the flux quantum and d is the interlayer dis- 
tance (~ 15 A in BSCCO). Indeed, in Fig. 7 the tri- 
angles demonstrate an upper bound for i? cl (ZD—Here, 
we take 2100 Aj-asj-a lower bound for \i("Pm an d 
10 fj,m for A c (0)E3lJ. We use the temperature depen- 
dence for AX ab (T) from o£ 6 (T) = l/[uj^\ 2 ab (T)] shown 
in the inset of Fig. 3a, and AA C (T) (squares in Fig. 2). It 
is clearly seen that neither the absolute value (too small 
compared to the experimental data) nor the temperature 
dependence (quasi- linear) agrees with the Hj{T) data. 

It is therefore quite natural to assume that a Bean- 
Livingston surface barrierEU impedes magnetic flux pen- 
etration into the sample and yields a higher entry field 
Hsb (T) . This assumption is also supported by the it- 
reversible behavior of the dissipation upon flux exited. 
In anisotropic superconductors in the quasi-2D regime, 
which holds in BSCCO up to temperatures very close to 
T c , the field Hsb(T) was shown to-be related only to the 
c-axis penetration length throughEj : 



H SB (T) 



47rA c (T)d' 



(6) 



The surface barrier might thus account for the ob- 
served value of the penetration field. Also, since Hg B (T) 
grows as 1/A C (T) (instead of l/A Q fcA c (T)), it is expected 
that the temperature dependence could show a better 



agreement. So we derive from the Hj{T) data an effec- 
tive penetration depth Aj(T) = <fio /[Air Hj(T)d] on the 
analogy of Eq. (6). The data are shown in Fig. 8 (full 
squares). We then try to determine A c (0) = Aj(T) — 
AA C (T) so as to fit \j(T) using the measured by ac- 
susceptibility and cavity perturbation technique AA C (T) 
(full circles in Fig. 8 extend the lower curve in Fig. 2 to 
higher temperatures). We find that both sets of data, 
namely Xj(T) and AA C (T), cannot be reconciled for any 
value we may assume for A c (0) in the entire temperature 
range. Therefore, the interpretation cannot be so simple. 
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FIG. 8. Plot of the temperature variation AA C (T) (full cir- 
cles, right scale) and of the effective length Aj(T) which is as- 
sociated to a surface barrier for the penetration of Josephson 
vortices (full squares, left scale) in sample #1. Open circles 
display the best fit using A c (0) = 35 /im and a scaling factor 
/3 = 6. Open symbols show the best fits using A c (0) = 20 /im 
(down triangles) and A c (0) = 100 /im (up triangles). 

More detailed quantitative investigation and calcula- 
tions of the penetration field in the presence of edge or 
surface defects, have led us to the conclusion that HjlT) 
is eventually controlled by such surface irregularitiesEil. 

We have carried out several checks (cleaving sample, 
cutting deep and wide grooves at the a6-plane, rotation 
crystal along it's c-axis) and established that large de- 
fects, greater or comparable with the penetration depth, 
within the a6-planes or at the edges of the crystal do not 
destroy the surface barrier (no significant change in the 
onset field of the microwave absorption was observed). 
However, if the dimensions of the surface irregularities 
are smaller than A then their presence can strongly in- 
fluence the screening current distribution. Indeed, the 
entrance field is deduced from the balance between the 
vortex attraction to the surface and the pushing force ex- 
erted by the screening Gwpent at the minimum distance 
£ (the vortex core size)EiH23. Near a small scratch the 



G 



current density can be many times larger than near the 
flat surface. This may substantially increase the force 
pushing vortices inside the superconductor and then de- 
crease the surface barrier and the entrance field. As 
shown in RefEU, in the specific case of a thin groove with 
the depth b > d at the surface, the penetration field 
Hj(T) = H SB (T)/P, where parameter f3 = (b/d) 1 / 2 > 1. 
Based on this estimation, we have attempted to fit our 
data on 
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AnH SB (T)d 



= /3[AA c (T)+A c (0)], (7) 



using two adjustable parameters f3 and A c (0) (see Fig. 8). 
We have obtained the best fit at [3 = 6 (or b ~ 500 A) 
and A c (0) = (35 ±15) /xm in Eq. (7), which is in agree- 
ment with our measurements by the first two methods. 
We also show in Fig. 8 smaller and larger values for A c (0). 
They allow us to set the uncertainty about our determi- 
nation of the penetration depth. We can also fit the data 
in the presence of a slit at the edge of the crystalo. The 
depth of the edge slit should be of the order of 10 /im, 
which is still small with respect to A c (0). The key result 
in this latter case is that it yields the same absolute value 
for A c (0). 



VI. CONCLUSION 

In conclusion, we have used three high-frequency tech- 
niques in studying anisotropic properties of BSCCO sin- 
gle crystals. The results obtained by different techniques 
are in reasonable agreement. We have observed almost 
linear dependences AA C (T) oc T at low temperatures. We 
have also determined the absolute value of A c (0), which 
is a factor of three higher in the optimally doped BSCCO 
sample than in the overdoped crystal. The ratio between 
the slopes of curves of AA C (T) in the range T <C T c is the 
same. These facts could be put down to dependences of 
A c (0) and AA C (T) on the oxygen content in these sam- 
ples. At the same time, the set of our experiments sug- 
gest very strongly the influence of defects in the samples, 
which relates the first penetration field Hj(T) of Joseph- 
son vortices to the c-axis penetration depth. In order 
to draw ultimate conclusions concerning the nature of 
the transport properties along the c-axis in BSCCO sin- 
gle crystals, studies of more samples with various oxygen 
contents are needed. 



VII. ACKNOWLEDGEMENTS 

This work was supported by the Centre National de 
la Recherche Scientifique - Russian Academy of Sciences 
cooperation program 4985, and by CREST and Grant- 
in- Aid for Scientific Research from the Ministry of Edu- 
cation, Science, Sports and Culture (Japan). The work 
at ISSP was also supported by the Russian Foundation 



1 M. R. Trunin, Physics- Uspekhi 41, 843 (1998); J. Super- 
conductivity 11, 381 (1998). 

2 T. Shibauchi, H. Kitano, K. Uchinokura, A. Maeda, 
T. Kimura, and K. Kishio, Phys. Rev. Lett. 72, 2263 (1994). 

3 J. Mao, D. H. Wu, J. L. Peng, R. L. Greene, and S. M. An- 
lage, Phys. Rev. B 51, 3316 (1995). 

4 H. Kitano, T. Shibauchi, K. Uchinokura, A. Maeda, 
H. Asaoka, and H. Takei, Phys. Rev. B 51, 1401 (1995). 

° D. A. Bonn, S. Kamal, K. Zhang, R. Liang, and 
W. N. Hardy, J. Phys. Chem. Solids 56, 1941 (1995). 

6 T. Jacobs, S. Sridhar, Q. Li, G. D. Gu, and N. Koshizuka, 
Phys. Rev. Lett. 75, 4516 (1995). 

7 T. Shibauchi, N. Katase, T. Tamegai, and K. Uchinokura, 
Physica C 264, 227 (1996). 

8 H. Srikanth, Z. Zhai, S. Sridhar, and A. Erb, J. Phys. 

Chem. Solids 59, 2105 (1998). 
9 H. Kitano, T. Hanaguri, and A. Maeda, Phys. Rev. B 57, 

10946 (1998). 

10 A. Hosseini, S. Kamal, DA. Bonn, R. Liang, and 
W. N. Hardy, Phys. Rev. Lett. 81, 1298 (1998). 

11 A. Buzdin and D. Feinberg, Phys. Lett. A165, 281 (1992). 

12 N. E. Hussey, A. Carrington, J. R. Cooper, and D. C. Sin- 
clair, Phys. Rev. B 50, 13073 (1994). 

13 L. N. Bulaevskii, V. L. Pokrovsky, and M. P. Maley, Phys. 
Rev. Lett. 76, 1719 (1996). 

A. E. Koshelev.P/M/s. Rev. Lett. 76, 1340 (1996). 

Matsuda, M. Gaifullin, K. Kumagai, M. Kosugi, and 
Hirata, Phys. Rev. Lett. 8, 1972 (1997). 
Sonin, Phys. Rev. Lett. 79, 3732 (1997). 
E. Gough and N. J. Exon, Phys. Rev. B 50, 488 (1994). 
V. Shovkun, M. R. Trunin, A. A. Zhukov, Yu. A. Nefy- 
odov, H. Enriquez, N. Bontemps, A. Buzdin, M. Daumens, 
and T. Tamegai, JETP. Lett. 71, 92, 2000. 
We note that 5f(T) includes the frequency shift due to the 
sample thermal expansion, which is essential for T > 0.7 T c 

in the T-orientationu. 

M. R. Trunin, Yu. A. Nefyodov, and H. J. Fink, Icond 



14 

15 y 

K 

16 E. 

17 C 

18 



D. 



mat/991121l| (submitted to JETP). 
H. Enriquez, N. Bontemps, A. A. Zhukov, D.V. 
M.R. Trunin, A. Buzdin, M. Daumens, and T. 



Shovkun, 
Tamegai, 
submitted. 

N. Bontemps, H. Enriquez, and Y. DeWilde, Physics 
and Materials Science of Vortex states, Flux Pinning and 
Dynamics, NATO Science Series 356 Ed. R. Kossovsky, 
S. Bose, W. Pan, and Z. Durusoy 387 (1999). 
' J. R. Clem, M. W. Coffey, and Z. Hao, Phys. Rev. B 44, 
2732 (1991). 

A. Schilling, F. Hulliger, and H. T. Ott, Z. Phys. B 82, 9 
(1991). 

S-F. Lee, D. C. Morgan, R. J. Ormeno, D. M. Broun, 
R. A. Doyle, and J. R. Waldram, Phys. Rev. Lett. 77, 735 



7 



(1996) . 

T. Tamegai, M. Sato, A. Mashio, T. Shibauchi, and S. 
Ooi, Proceedings of the 9th International Symposium on 
Superconductivity, Advances in Superconductivity IX, Ed. 
S. Nakajima and M. Murakami, 621, Springer- Verlag Tokyo 

(1997) . 

N. Nakamura, G. D. Gu, and N. Koshizuka, Phys. Rev. 
Lett. 71, 915 (1993). 

F. Steinmeyer, R. Kleiner, P. Miiller, H. Miiller, and 



K. Winzer, Europhys. Lett. 25, 459 (1994). 

O.K.C. Tsui, N. P. Ong, and J. B. Peterson, Phys. Rev. 

Lett. 76, 819 (1996). 

B. Khaykovich, E. Zeldov, D. Majer, T. W. Li, P. Kes, and 
M. Konczykowski, Phys. Rev. Lett. 76, 2555 (1996). 

C. P. Bean and J. B. Livingston, Phys. Rev. Lett. 12, 14 
(1964). 

P. G. De Gennes, Superconductivity of Metals and Alloys, 
Benjamin, NY, (1996). 



8 



